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Effects of single-site anisotropy on mixed diamond chains with spins 1 and 1/2 are investi- 
gated in the ground states and at finite temperatures. There are phases where the ground state 
is a spin cluster solid, i.e., an array of uncorrelated spin-1 clusters separated by singlet dimers. 
The ground state is nonmagnetic for the easy-plane anisotropy, while it is paramagnetic for the 
easy-axis anisotropy. Also, there are the Neel, Haldane, and large- Z) phases, where the ground 
state is a single spin cluster of infinite size and the system is equivalent to the spin-1 Heisenberg 
chain with alternating anisotropy. The longitudinal and transverse susceptibilities and entropy 
are calculated at finite temperatures in the spin-cluster-solid phases. Their low-temperature 
behaviors are sensitive to anisotropy. 

KEYWORDS: mixed diamond chain, anisotropy, frustration, spin cluster solid, Neel phase, Haldane phase, 
large- D phase 



1. Introduction 

Frustration in quantum spin systems has been one of 
the most exciting subjects studied in the field of mag- 
netism over these past decades. 1 ' 2 ) A number of theo- 
retical and experimental investigations have revealed ex- 
otic phenomena induced by the interplay of frustration 
and quantum fluctuation. In the theoretical approach, 
many exactly solvable models have played crucial roles 
in elucidating frustrated quantum magnetism. As typi- 
cal examples, there exists a class of models whose ground 
states are exact spin cluster solid (SCS) states because of 
frustration. The SCS state is defined as a tensor product 
state of exact local eigenstates of cluster spins. A dimer 
state is the simplest type of SCS state. For example, 
the Majumdar- Ghosh model 3 ) has a dimer ground state, 
which is a prototype of spontaneously dimerized states 
in onc-dimcnsional frustrated mag nets. 4) The Shastry- 
Sutherland model, 5 ) which corresponds to the material 
SrCu 2 (B03)2, 6 ' 7 ) also has a dimer ground state. 

The diamond chain is a frustrated spin chain with 
exact SCS ground states that are different from dimer 
states. The lattice structure of the diamond chain is 
shown in Fig. 1. In a unit cell, there are two kinds of 
nonequivalent lattice sites occupied by spins with mag- 
nitudes S and t; we denote the set of magnitudes by (S, 
t). One of the authors and coworkers 8 ' 9 ) introduced this 
lattice structure and generally investigated the case of 
(S, S), i.e., the pure diamond chain (PDC). After that, 
Niggcmann et al. 10 ' 11 ' argued about a series of diamond 
chains with (S, 1/2). 

The mixed diamond chain (MDC) is defined as a dia- 
mond chain with (S, 5/2) for the integer 5. 12 ) Recently, 
extensive investigation on the MDC has been carried out 
by the present authors. 12 ~ 15 ) The MDC is of special in- 
terest among diamond chains, because only the MDC has 
the Haldane phase in the absence of frustration so that 
we can observe the transition from the Haldane phase 
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Fig. 1. Structure of the diamond chain. Spin magnitudes in a unit 
cell are indicated by S and t; we denote the set of magnitudes 
by (S, t). The PDC is the case of S = r, while the MDC is the 
case of S = 2t with integer or half-odd integer r. 



to a SCS phase induced by frustration. In contrast, di- 
amond chains of other types have ferrimagnetic ground 
states for weak frustration. 

The common features of any types of diamond chains 
are that they have an infinite number of local conser- 
vation laws and more than two different types of ex- 
act SCS ground states are realized depending on the 
strength of frustration. The MDC with (1,1/2) has 3 dif- 
ferent paramagnetic SCS phases accompanied by spon- 
taneous translational symmetry breakdown (STSB) and 
one paramagnetic phase without STSB in addition to a 
nonmagnetic Haldane ground state in a less frustrated 
region. 10,12 ) The SCS structures of the ground states are 
also reflected in characteristic thermal properties, as re- 
ported in ref. 13. 

Thus, the diamond chain serves as a playground for in- 
vestigating various additional effects on frustrated mag- 
netism on a well-founded basis. Therefore, various mod- 
ifications of the PDC and MDC have been examined by 
many researchers. 14 ~ 19 ) It is found that the natural min- 
eral azurite consists of distorted PDCs with spin 1/2 and 
the magnetic properties of this material have been exper- 
imentally studied in detail. 20 ' 21 ) Other materials have 
also been reported. 22 ' 23 ) In addition, as reviewed in ref. 
14, the MDC is related to many other important models 
of frustrated magnetism. 24 ~ 33 ) 

In the present work, we concentrate on the effects 
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of the single-site anisotropy D on the ground-state and 
finite-temperature properties of the MDC with (1,1/2). 
In many quantum spin systems, it is known that the D- 
term changes the ground state drastically. A well-known 
example is the effect of the Z?-term on the Haldane state 
of spin-1 antifcrromagnetic Heisenberg chains. A large 
positive D term destroys the hidden order in the Hal- 
dane phase and drives the ground state into the large- D 
phase that has no specific order. On the other hand, a 
large negative £>-term stabilizes the Neel order. 34 ~ 39 ' The 
D-terms in the higher spin and mixed spin Heisenberg 
chains induce rich phase diagrams. 39-42 ) The alternating 
D-terms in spin-1 Heisenberg chains can pin the fluctu- 
ating hidden order to induce a long-period Noel order. 43 ' 
They can also induce the partial ferrimagnctic order in 
higher-spin chains. 44 ' Considering the high degeneracy 
of the ground states of the MDC, 12 ' we may expect an 
even richer phase diagram in the MDC with the single- 
site anisotropy. 

Thus far, despite the theoretical relevance of the MDC, 
no materials described by the MDC have been found. 
Nevertheless, synthesizing MDC materials is not an un- 
realistic expectation in view of the successful synthesis of 
many novel magnetic materials such as molecule-based 
magnetic materials. 45 ' In general, it is natural to ex- 
pect the single-site anisotropy on the spin-1 site in real 
MDC compounds. From this viewpoint, it is important 
to present theoretical predictions on the ground state 
and finite-temperature properties of anisotropic MDCs 
to widen the range of candidate materials of the MDC 
and raise the possibility of their synthesis. We find that 
the effect of the single-site anisotropy not only gives 
quantitative correction to physical quantities but also 
changes the ground state and low-temperature behavior 
qualitatively. Among them, the SCS ground states are 
found to be sensitive to such anisotropy. 

This paper is organized as follows. In §2, the Hamilto- 
nian for the MDC with the single-site anisotropy is pre- 
sented. In §3, the MDC with the single-site anisotropy 
is examined by numerical methods and the ground-state 
phase diagram is obtained. Various limiting cases in the 
phase diagram are discussed by the perturbation method. 
In §4, the finite-temperature behaviors of the longitudi- 
nal and transverse susceptibilities and entropy are in- 
vestigated numerically and analytically by extending the 
method developed in ref. 13. The last section is devoted 
to the summary and discussion. 

2. Hamiltonian 

We consider the MDC with the single-site anisotropy 
on spin-1 sites described by the Hamiltonian 



H = 



i=i 



Si + S 



i+i) 



(r 



(i) 



.(2), 



(1)^(2) 



(1) 



1=1 



where Si is the spin-1 operator, and r ; and t| 2 ^ are the 
spin-1/2 operators in the ith unit cell. The total number 
of unit cells is denoted by L. In the case of D = 0, eq. (1) 



reduces to the Hamiltonian of the isotropic MDC. 12 ' In 
what follows, we set the energy unit as J = 1. 

Before proceeding, we analyze the classical version of 
the quantum Hamiltonian (1). We obtain the classical 
ground state for any A and D in Appendix. We have four 
classical phases separated by the phase boundaries A = 2 
and D = 0. In the two phases of A > 2, the ground-state 
spin configurations can be locally modified without an 
energy increase. This classical situation corresponds to 
the quantum situation that there are an infinite number 
of low-energy states that arc transformed from each other 
by local modification. Such quasi-degenerate low-energy 
states may enhance quantum fluctuations to contribute 
to the appearance of exotic quantum states for large A 
in the quantum system even for finite D. 

3. Ground-State Phase Diagram 

3. 1 General formulation for phase boundaries 

The Hamiltonian (1) has a series of conservation laws. 
To see it, we rewrite eq. (1) in the form 



1=1 



(S l T l +T l S l+1 ) + 



A 



i=i 



where the composite spin operators T; are defined as 



+ r[ 2) (Z = 1,2,-..L). 



Then, it is evident that 

[T 2 ,^]=0 (Z = l,2,---Z). 



(2) 

s 
(3) 

(4) 



Thus, we have L conserved quantities T 2 for all I, even 
for D ^ 0. By defining the magnitude T; of Tj by 
T\ = Ti(Ti + 1), we have a set of good quantum numbers 
{Ti; I = 1, 2, L}. Each T; takes a value of or 1. The 
total Hilbert space of the Hamiltonian (2) consists of sep- 
arate subspaces, each of which is specified by a definite 
set of {Ti}, i.e., a sequence of and 1. 

A spin pair of TJ = is a singlet dimer, that cuts off 
the correlation between Si's at both sides as seen from 
eq. (2). Hence, when a segment is bounded by two X} = 
pairs, it is isolated from other parts of the spin chain. 
The segment including n successive Ti 's with XJ = 1 and 
n+ 1 Si's is called a cluster- n as in the isotropic case. 12 ' 
A cluster-n is equivalent to an antifcrromagnetic Heisen- 
berg chain consisting of In + 1 effective spins with mag- 
nitude 1 with alternating anisotropy. The Hamiltonian is 
written as 

2n n+1 

1=1 1=1 

where S 2 i-i = Si and S 2 i = T/. 

The ground state of the total diamond chain is an- 
alyzed by describing each cluster by the Heisenberg 
model (5). For D = 0, we have obtained the complete 
ground-state phase diagram in ref. 12. In particular, 
we found successive phase transitions between dimer- 
cluster-n (DCn) phases with n = 0, 1, 2, 3, and oo as 
A decreases. The DC0 phase is also called the dimer 
monomer phase where 7] = for all I. The DCn phase 
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with n = 1, 2, or 3 is the phase where the ground state 
consists of an alternating array of cluster-n's and sin- 
glet dimers. The DCoo phase is the Haldane phase for 
the equivalent Hciscnbcrg chain with T; = 1 for all I. 
The description in terms of cluster-n and DCn states is 
naturally extended to the finite- D regimes. By the same 
argument as that in ref. 12, the phase boundary between 
the DC(n — 1) and DCn phases for finite D is given by 

A c (n - l,n;D) 

= (n+ l)E G {2n — 1,D) — nE G (2n + 1, D), (6) 

where Eq(2u + 1,D) is the ground-state energy of the 
Hamiltonian (5). These phase transitions are of the first 
order, since they take place as level crossings between two 
eigenstates of the original Hamiltonian (1) characterized 
by different sets of quantum numbers {T{\. The direct 
transition from the DCn phase to the DCoo phase takes 
place at 

A c (n,oo;L>) = E G (2n + l,D) - 2(n + l)e G (oo,£>), (7) 

if A c (n, oo; D)>A c (n, n + 1;D), where €q(oo 7 D) is the 
ground-state energy of the Hamiltonian (5) with n — > oo 
per unit cell. 

The above arguments are almost parallel to those in 
the isotropic case. For the physical conclusion, however, 
we observe a significant difference from the isotropic case, 
as described in the following sections. 

3.2 Ground- state phase diagram determined by numer- 
ical method 

To obtain the ground-state phase diagram using eqs. 
(6) and (7), we have numerically diagonalized the Hamil- 
tonian (5) to calculate the ground-state energy Ec(2n + 
1,-D). The ground-state energy per unit cell eg(cg,D) 
is calculated by the infinite-size DMRG by estimating 
the energy increment accompanied by an addition of a 
unit cell at the center of the chain. The phase diagram 
is shown in Fig. 2. Because of the presence of fine struc- 
tures of the phase diagram, we present a series of phase 
diagrams magnified appropriately. In the top figure, the 
broken and dotted lines are the results of the perturba- 
tion analysis described in the next section. 

The main features characteristic of the anisotropic case 
are as follows. 

(i) The total spin S to t of cluster-n is not a good quan- 
tum number, while its z-component remains a good 
quantum number. The ground state of cluster-n has 
S" t 2 ot = for D > and S t z ot = ±1 for D < 0. There- 
fore, the DCn phase with finite n is nonmagnetic for 
D > and paramagnetic for D < 0. 

(ii) The DCoo state corresponds to the ground state of 
a spin-1 Heisenberg chain that has the single-site 
anisotropy and D altcrnatingly. It is known that 
the ground-state phase diagram of this model con- 
sists of a Neel phase for D < —0.695, a Haldane 
phase for —0.695 < D < 3.28, and a largc-D phase 
for D > 3.28. 43) 

(iii) There appear DCn phases with n = 0, 1, 2, and 3 
if A is above the critical value, which depends on D. 
These phases continue to the isotropic DCn phases. 
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Fig. 2. Ground-state phase diagram of the anisotropic mixed di- 
amond chain. The figures are magnified appropriately from (a) 
to (d). In (a) and (b), some phases are not indicated because 
they are invisibly narrow on the scale of these figures. In (a), 
the analytical results in §3.3 are also shown by the dotted and 
broken lines. 
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(iv) There appears another DC3 phase separated from 
the line of D = and in contact with the end-point 
of the Neel-Haldanc transition line in the DCoo re- 
gion. 

(v) The critical values of A between the DCn phases 
decrease with \D\ from those in the isotropic case. 

(vi) For sufficiently large values of \D\ , only the DCO and 
DCoo phases exist. 

3.3 Limiting cases 

3.3.1 DC0-DC1 phase boundary 

The phase boundary between the DCO and DC1 phases 
is obtained by substituting the eigenvalues E G (1, D) and 
-Eg (3, D) into eq. (6). These eigenvalues can be analyt- 
ically determined. Obviously, E G (1,D) = for D > 
and E G (1, D) = D for D < 0. 

For D > 0, the ground state of a cluster-n belongs 



to the subspace of S*, 



0, P = 1, and T = 1, where 



S* ot is the z-componcnt of the total spin. P is the space 
inversion parity, and T is the time reversal parity. Within 
this subspace, E G (3,D) is the smallest solution of the 
following eigenvalue equation: 

E G (3,D) 3 -(3D -1)E G (3, D) 2 



2(D 2 -D- 3)J5 G (3, D) + 8D = 0. 



(8) 



Therefore, A c (0, 1) = -E G (3, D). For small values of D, 
this implies A c (0, 1) ~ 3 - ffD. 

For D < 0, the ground state of a cluster-n belongs to 
the subspace of S* ot = ±1 and P = 1. Within this sub- 
space, E G (3,D) is the smallest solution of the following 
eigenvalue equation: 

E G (3, D) A — (5D — 2)E G (3, D) 3 

+ (SD 2 -6D -b)E G (3,D) 2 

+ (-AD 3 + 4D 2 + IAD - 6)E G (3, D) 

- 8D(D - 1) = 0. (9) 

Therefore, A c (0, 1) = 2D - E G (3, D). For small values of 
\D\, this implies A c (0, 1) ~ 3 + 

3.3.2 D -¥ -oo 

For D < and \D\ 1, we apply the degenerate per- 
turbation theory to the ground-state energy of a cluster- 
n. We regard the diagonal part of H n as the unperturbed 
Hamiltonian and the off-diagonal part as the per- 
turbation Hamiltonian Hn as follows: 

n+l n 



1=1 



i=i 



^ 1) = f E[^ ++5 m) T r+ h - c -] 



(10) 



1=1 



The expansion parameter a is introduced, which will be 
set equal to unity after the calculation. Up to the third- 
order perturbation calculation in a, we have 



E G (2n + !,£>)- - (n + 1)|D| - 2n 



(2n - 2) 

W\ 



-or 



2 2 (3n-3) o 



\D\ + 2 2(|D|+3) 2 [\D\+2f 



(11) 



for n > 1. Substituting this expression into eq. (6), the 
phase boundaries are determined up to the third order 
in a as follows. 

(i) The DCn-DC(n + 1) transition point for n > 1 is 
given by 



A c (n, n + 1; D) ~ 2 



4 2 2 2 
-a — ra 



+ 



|£>| + 3 |D|+2 
1 



y\D\+3y~ (\D\+2)i a3 - (12) 

The n independence of eq. (12) means the direct 
transition between the DC1 and DCoo phases within 
this approximation. 

(ii) The DC0-DC1 transition point is given by 



A c (0, !;£>)- 2 1 



1 



1 



\D\+2 2(|D| + 2) 



(13) 



where the exact relation -Eb(l) D) = —\D\ is used. 

Setting a = 1, these approximate phase boundaries are 
drawn by the dotted lines in Fig. 2(a). 

3.3.3 D -> oo 

In this limit, all the 5-spins are in the state S z = 0. 
Therefore, we find the following: 

3L 

(i) The ground-state energy of the DCO phase is — —A. 

(ii) In the DCoo phase, the effective interaction between 
T-spins is 0(1/ D) as described in ref. 43. The effec- 
tive Hamiltonian is given by 

H e S = ~tA 



E 

i=i 



eff 



D, 



[J eS (TfT? +1 + T?T? +l ) + D cS (T? 2 - 2)] , 

(14) 
(15) 



err 



D 



The ground-state energy of this effective model is 
calculated by the infinite-size DMRG method as 

— — 2.6995J G ff per site. According to the numeri- 
cal calculation, there appear no DCn phases with 
1 < n < oo for D > 1.24. Taking this into account, 
we conclude that the direct transition between DCO 
and DCoo phases takes place at 

A c (0,oo ;j D)^^ (16) 

for sufficiently large values of D up to 0(l/D). This 
approximate phase boundary is drawn by the broken 
line in Fig. 2(a). 
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4. Finite- Temperature Properties 

4-1 General formula 

Following the method described in ref. 13, we can cal- 
culate the thermodynamic expectation value of the ex- 
tensive quantity Q per unit cell as 

(Q) _ (Qci) 

L (Ld) ' 

where Q c \ is the physical quantity Q for each cluster-n 
and L c \ = n + 1 is the number of unit cells per cluster- 
n. The average (Q c i) of Q c \ over the grand canonical 
ensemble is written as 



(17) 



2J Z c i(n) (Q c \(n)) c 



D /3( M -A/4)(n+l) 



(<5ci> 



(18) 



B /3( M -A/4)(n+l) 



n=0 



where n is the chemical potential of a cluster-n. 
(<3ci("-)) can and Z c i(n) are the canonical average of the 
physical quantity Q and the partition function of a 
cluster-?i with fixed n, respectively. These are defined 

by 

(Qc 1 N) can = ^-EE^ (2 " +1 ^ l) 

x (2n + 1, v\ M5| Qci |2n + 1, i/; AQ) , (19) 

(20) 



W = ^E 



-0B(2n+l,iy;Af c i: 1 ) 

' 1 



where i5(2n + 1, M*j) and |2n + 1, i/; M^) are, respec- 
tively, the f-th eigenenergy and eigenstate of the spin-1 
chain with length 2n + 1 and magnetization M^. The 
chemical potential u is determined by the condition 



^e^Z d (n)e*- A / 4 »("+ l )=l. 

n=0 



(21) 



^.J? Formulae for entropy and magnetic susceptibility 

The entropy per unit cell is calculated using the same 
formula as that used in the isotropic case derived in ref. 
13 as 



(22) 



where E c \ = E + (n — 3)A/4 is the energy per cluster-n. 

The longitudinal magnetic susceptibility x\\ is calcu- 
lated using the same formula as that used in the isotropic 
case by the direct application of eqs. (17), (18), and (20). 
This gives 



All 



1 {Ml) 



(23) 



T (La) ' 

The transverse susceptibility x± 1S calculated as dis- 
cussed below. 

The contribution to the Hamiltonian due to the trans- 
verse field is given by 

L 

U x = -Hj2(Sf + Tn- (24) 
l=i 



Because each T ; commutes with H x , all Ti's are con- 
served even in the presence of the transverse field. There- 
fore, the magnetic susceptibility of an anisotropic MDC 
(1) is the sum of the contributions of cluster-n's. 

Using the standard linear response theory, the canoni- 
cal expectation value of the transverse magnetic suscep- 
tibility of a cluster-n is given by 

" v ' J1P 1= 0i/i,i/ 2 
e -PE{2n+l,UT_;M'\) _ e -0E(2n+l,u 2 ;M^+l) 

X E(2n + 1, v 2 - M\ + 1) - E(2n + 1, v V ,M^) 

x |(2n+l,z/ 2 ;M5 + l|M c ||2n+l, i y 1 ;M5)| 2 , (25) 

where E(2n + l,u;M^) and |2n + l,u;M^) are the v- 
th eigenenergy and eigenstate of the spin-1 chain with 
length 2n + 1 and magnetization M*j. Thus, the trans- 
verse susceptibility of the MDC per unit cell is given by 

(Xcix) 



X± = 



(Lei) ' 



(26) 



where (• ■ ■) means the grand canonical average defined 
in eq. (18). 

In the DCn ground state, this formula reduces to 

1 1 

X± (n + 1) E(2n + 1, i/; 1) - E G (2n + 1,0) 

x |<2n+l,i/;l|A/+|2n+l,G;0)| 2 (27) 
for D > and to 

I 

2(n + 1) ^ I FA 



X± 



-E(2n+l,u;2) - E G (2n+l,l) 

x |(2n + l,^;2|A/+|2n+l,G;l)| 2 
1 

E(2n + 1, u; 0) - E G (2n + 1, -1) 

x |(2n + l,^;0|M+|2n+l,G;-l)| 2 ] (28) 

for D < 0. Here, 2n + 1, G; M c \) is the ground state 
of a cluster-n with magnetization M c \. For D < 0, the 
ground states are degenerate, so that the summation is 
taken over all the ground states. It should be noted that 
Xi_ remains finite at T = in both cases. 

4-3 Results for entropy and magnetic susceptibility 

In cq. (18), the summation over n is taken over all non- 
negative integers. In the actual numerical calculation of 
finite-temperature properties, however, we can only in- 
clude the contribution of cluster-n's with finite n. Hence, 
we cannot expect reliable results if the ground state is 
the DCoo state. Therefore, we limit ourselves to the pa- 
rameter region with DCn ground states with finite n. 
In what follows, we present the results for the magnetic 
susceptibility and entropy calculated including the con- 
tribution of cluster-n's with < n < 5. The error due to 
this cutoff procedure is estimated for the isotropic case 
in ref. 13. Because this estimation is based on an entropic 
argument in the high-temperature region, it is also valid 



6 J. Phys. Soc. Jpn. 



Full Paper 



K. Hida and K. Takano 



in the present case. Thus, we expect that the missing 
entropy is within 3% of the total entropy. 
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Fig. 3. Temperature dependences of (a) longitudinal magnetic 
susceptibility X\\ times T, (b) transverse magnetic susceptibil- 
ity x±> an d ( c ) entropy S for D = —0.01. The chosen values of 
A are in the DCra phases and also at the DC(n + l)-DCn phase 
boundaries. The large filled right triangles on the ordinates are 
the ground-state values in each phase and the large open ones are 
those at the phase boundaries. The insets are magnified figures 
in the low-temperature regimes. 



To demonstrate the sensitivity of the low-temperature 
behavior to the anisotropy D, we present the numerical 
results for D = ±0.01 as representatives of the easy-plane 
and easy-axis anisotropics. For a negative £>, the longi- 
tudinal and transverse susceptibilities are shown in Figs. 
3(a) and 3(b), respectively, and the entropy is shown 
in Fig. 3(c). For a positive D, both the susceptibilities 
are shown in Fig. 4(a) and the entropy is shown in Fig. 
4(b). The difference between the low-temperature behav- 
iors for the positive and negative D's is distinct even 



for a small \D\ = 0.01. In the low-temperature limit, 
the longitudinal magnetic susceptibility xil an d the en- 
tropy S tend to zero for D > 0, except at the ground- 
state phase boundary. On the other hand, they behave 

as X\\T — and S ~ for D < if the ground 

11 n + 1 n + 1 

state belongs to the DCn phase. The low-temperature 

limiting values of x± are calculated using eqs. (27) and 

(28). These limiting values are shown by the large filled 

right triangles on the ordinates in Figs. 3 and 4. 

At the DC(n— l)-DCn phase boundary, cluster-(n — 1) 

and cluster-n coexist. In this case, the residual entropy 

can be estimated by the combinatory argument similarly 

to that in the isotropic case. If we denote the number 

fraction of cluster- (n — l)'s by x and that of cluster-n's 

by 1 — x, then the entropy can be estimated as 



S = -- 



L 



[xlax + (1 — x) ln(l — x) — In g] , 
it ~r x — ju 

(29) 

where g is the degeneracy of the ground state of cluster- 
n, i.e., g = 1 for D > 0, g = 2 for D < 0, and g = 3 for 
D = 0. By optimizing S with respect to x, we find 



.9 



1 -x 



(30) 



Using x that satisfies eq. (30), the entropy and suscepti- 
bility can be estimated as 



S , x 
— = In , 

L 1-x' 

X\\,X xXd\\,A n -!) + (!- x )Xd\\,A n ) 



(31) 



(32) 



L nx+ (n + 1)(1 -x) 

It should be noted that the residual entropy remains fi- 
nite at the phase boundary even for D > owing to 
the mixing entropy of cluster- (n — 1) and cluster-n. For 
D < 0, Xci || = 1/r f° r ah n - Therefore, we have 

(33) 



L 



T(n + 1 



For D > 0. 



Xll 

— = 0. To calculate x±i we have to esti- 
L 

mate Xci\\,±i n ~ 1) anc l Xci\\,±i n ) numerically using eqs. 
(27) and (28). These low-temperature limiting values are 
shown by large open right triangles on the ordinates in 
Figs. 3 and 4. 

5. Summary and Discussion 

We examined the anisotropic mixed diamond chain 
with spins 1 and 1/2 that has the single-site anisotropy 
D on spin-1 sites. In the ground-state phase diagram 
(Fig. 2), there are DCn phases with n = 0, 1, 2, and 
3 if A is above the critical value, which depends on D. 
These phases continue to the isotropic DCn phases ex- 
cept for a DC3 phase separated from the line of D = 0. 
As in the isotropic case, the DCn ground state with fi- 
nite n is an alternating array of finite-length cluster-n's 
and singlet dimers. It is nonmagnetic for the easy-plane 
anisotropy (Z?>0), and is paramagnetic for the easy-axis 
anisotropy (D<0). For smaller A, the ground-state phase 
is one of the three DCco phases: i.e., the Neel phase for 
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Fig. 4. Temperature dependences of (a) longitudinal (x||) an d 
transverse (x±) magnetic susceptibilities and (b) entropy S for 
D = 0.01 . The chosen values of A are in the DCn phases and also 
at DC(n + l)-DCn phase boundaries. The large filled right trian- 
gles on the ordinates are the ground-state values in each phase 
and the large open ones are those at the phase boundaries. 



D < -0.695, the Haldane phase for -0.695 < D < 3.28, 
and the large- D phase for D > 3.28. The critical val- 
ues of A between the DCn phases decrease with \D\ from 
those in the isotropic case. For sufficiently large values 
of \D\, only the DCO and DCoo phases exist. The phase 
boundary between them tends to A =2 for D —> — oo and 
to A =0 for D -> oo. 

The above features of the ground-state phase diagram 
are qualitatively understood, if we incorporate the effect 
of quantum fluctuation into the classical phase diagram 
explained in Appendix. In the two classical phases with 



A > 2, a pair of Tj and can locally rotate around 
the fixed axis determined by Si's without increasing the 
energy, as seen in Appendix. The anisotropic term re- 
stricts the possible direction of the fixed axis: for D < 0, 
the fixed axis is only parallel or antiparallel to the z-axis, 
and for D > it is confined in the xy plane. The quan- 
tum effect of D corresponding to this classical restriction 
is to reduce the quantum fluctuation of Si 's only. There- 
fore, for A > 2, the anisotropic term does not suppress 
the quantum fluctuation corresponding to the classical 
local rotations of t| 's around the fixed axis of Si . On 
the other hand, for A < 2, the classical ground state is 
antiferromagnetically ordered so that the quantum fluc- 
tuation is collectively concerned with all the spins and is 



(2) 



fore, the D-term suppresses the collective quantum fluc- 
tuation, and hinders the energy gain by the quantum 
fluctuation. This explains the reason why the region of 
the exotic ground state for large A extends to lower values 
of A with an increase in the magnitude of the anisotropy 
\D\. 

The features of the ground-state phases and phase 
transitions are reflected in physical quantities at low 
temperatures. The temperature dependences of entropy, 
and longitudinal and transverse susceptibilities are cal- 
culated. The low-temperature behavior is sensitive to the 
sign of D for large A, reflecting the change of the ground- 
state phase. The transverse susceptibility remains finite 
even at T = 0, as long as D ^ 0. 

Generally, the introduction of anisotropy to a phys- 
ical model increases the possibility that a correspond- 
ing material is realized. If an anisotropic MDC mate- 
rial is synthesized, the anisotropy may be controlled by 
changing the crystal field that reflects the environment 
of spin-1 magnetic ions. This would be accomplished by, 
e.g., applying pressure and/or changing the nonmagnetic 
ions around the magnetic ones. In contrast to the spin- 
1 chain with the Haldane gap, the ground state of the 
MDC is sensitive to the weak anisotropy. Therefore, the 
paramagnetic-nonmagnetic transition should take place 
with a small variation in the crystal field. It would also 
be possible to observe the quantum phase transitions be- 
tween the DCn states with different spatial periodicities. 

In refs. 14 and 15, we have reported the effect of dis- 
tortion on the isotropic MDC. It is predicted that various 
types of distortion can induce various types of quantum 
phases such as the Haldane phase with STSB, and the 
quantized and partial ferrimagnctic phases. The inter- 
play of anisotropy and distortion would lead to an even 
wider variety of phenomena. Their investigation is left 
for future studies. 

The numerical diagonalization program is based on the 
package TITPACK ver.2 coded by H. Nishimori. The nu- 
merical computation in this work has been carried out 
using the facilities of the Supercomputer Center, Insti- 
tute for Solid State Physics, University of Tokyo, the Su- 
percomputing Division, Information Technology Center, 
University of Tokyo, and the Yukawa Institute Computer 
Facility, Kyoto University. KH is supported by a Grant- 
in-Aid for Scientific Research (C) (21540379) from the 
Japan Society for the Promotion of Science. 



Appendix 

We examine the phase diagram of the classical ver- 
sion of the quantum Hamiltonian (2), where Si and Tj 
are interpreted as classical vectors with length S and 
t = S/2, respectively. The classical Hamiltonian is then 
expressed in the following two forms: 



u 



clas 



2Ti + Si) 2 - 2(2 - A)T? - S, 



■m 



clas 



not restricted to individual pairs of r 



(i. 



and t\ 2 \ There- 



T l + Si) 2 -2Si 



clas 



(34) 
(35) 
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.(1) 



r| 2) , and Uf™ 



where Si = Si + S l+1 , T t 
Eti(DSf-XS 2 /4). 

In the case of A < 2, the expression (34) shows that 



W 



is minimized if \Si 



2S, IT, 



S, and \Ti 

(<*)> 



^Si\ = irrespective of D. Hence, all the Si's (t; 's) 
in the chain are aligned parallel (antiparallel) to a fixed 
axis, and the ground state is antiferromagnetic. The fixed 
axis is the z-axis for D < 0, while it may be any direction 
in the xy plane for D > 0. The ground state is elastic 
in both cases, since any local modification of the spin 
configuration increases the energy. 

In the case of A > 2, the expression (35) reveals that 
n clas is minimized if \Si\ = 25 and |Tj + S t /X\ = 
irrespective of D. Hence, Si's in the chain are aligned 
parallel to a fixed axis, and t, and r[ 2 ^ form a trian- 
gle with Si/X. The fixed axis is the z-axis for D < 0, 
while it may be any direction in the xy plane for D > 0. 
Then, the ground state has a local arbitrariness, because 
Tj and Tj may be rotated about the axis of Si and 
Si-\-\ without raising the energy. All the ground states 
are ferrimagnetic with a magnetization (1 — 2/X)SN. 

Thus, we have the classical phase diagram consisting 
of the four phases separated by the phase boundaries of 
A = 2 and D = 0. 
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